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■ ABSTRACT 

In this paper we describe the photometric cahbration of data taken with the near- 
' (— i '' infrared Wide Field Camera (WFCAM) on the United Kingdom Infrared Telescope 

p ,! (UKIRT). The broadband ZYJHK data are directly calibrated from 2MASS point 

I ■ sources which are abundant in every WFCAM pointing. We perform an analysis of 

Q I spatial systematics in the photometric calibration, both inter- and intra-detector and 

> . show that these are present at up to the ~5 per cent level in WFCAM. Although the 

^ ' causes of these systematics are not yet fully understood, a method for their removal 

is developed and tested. Following application of the correction procedure the photo- 
metric calibration of WFCAM is found to be accurate to ~ 1.5 per cent for the JHK 
' bands and 2 per cent for the ZY bands, meeting the survey requirements. We investi- 

gate the transformations between the 2MASS and WFCAM systems and find that the 
Z and Y calibration is sensitive to the effects of interstellar reddening for large values 
of E{B — Vy , but that the JHK filters remain largely unaffected. We measure a small 
correction to the WFCAM F-band photometry required to place WFCAM on a Vega 
system, and investigate WFCAM measurements of published standard stars from the 
list of UKIRT faint standards. Finally we present empirically determined throughput 

. measurements for WFCAM. 
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1 INTRODUCTION be uniform across the whole sky to better than 2 per cent 



This paper describes the photometric calibration of the 



accuracy (Nikolaev et al. 2000). 

/TTT^T-i^mx , n , , ^fals papeT we investigate the accuracy and homo- 
United Kmgdom Infrared Telescope (UKIRT) Wide Field ,,rn/^ATVT i-u i- u i • r j ^ i i 
„ ,x , , eeneity ot the Wl< CAM calibration, by analysis ot data taken 
Camera (WFCAM) ZYJHK passbands spanning 0.84 - n ^ ^ r ■ \- rj., 

^ , , r m the first two years of science operations, the survey re- 

2.37 um. The WFCAM photometric system, based on syn- . ^ ■ ^ • i i ^ ^ • i-i j_- • i i 

I , quirement is to provide a photometric calibration m broad 

thetic colours, is described in Hewett et al. (2006). , , r,,. ^ ^ r, ^ -xi i r i 

^ ' band filters accurate to 2 per cent with a goal of 1 per cent. 

WFCAM is currently mounted on UKIRT about 75 per rn ^ ^ i ii ii i-i. ^ j i, i ^ ii, 

■' ^ lo test whether the calibrated photometry passes the re- 
cent of the time, the majority of which is used to per- . , -j ^i, r n • ^ ^ 

' quirements, we consider the following tests: 
form a set of five public surveys, under the umbrella of the 

UKIRT Infrared Deep Sky Survey: UKIDSS, (Lawrence et • for repeat observations of non-variable stars (of sufli- 

al. 2007, Warren et al. 2007), as well as a number of cam- cient signal-to- noise) , the measured rms (root mean square) 

paign projects and other smaller programmes. should be no worse than 0.02 magnitudes. 

All WFCAM data are processed and calibrated by the • the calibration should be robust against the effects of 

Cambridge Astronomical Survey Unit (CASU) using an au- reddening, i.e. the calibrated frame zeropoints, measured on 

tomated pipeline (Irwin et al. 2008). Due to the large data photometric nights, should have a standard deviation < 0.02 

rate of WFCAM (150-230 Gigabytes of raw image data are magnitudes as a function of reddening. 

taken per night), a reliable and robust method for photomet- • the photometric calibration should be uniform across 

ric calibration is required. The primary photometric calibra- the sky, and not depend on the observed stellar population or 

tors are drawn from stars in the Two Micron all Sky Survey the number of stars in the image. Specifically, the calibrated 

(2MASS, Skrutskie et al. 2006), present in large numbers in frame zeropoints (on photometric nights) should have a la 

every WFCAM exposure (the median number is around 200 scatter of < 0.02 magnitudes even in sparse regions of sky, 

per pointing). The 2MASS calibration has been shown to e.g. at high Galactic latitude. 
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Table 1. Summary of photometric calibration corrections in- 
cluded in each UKIDSS release from the WSA. The columns 
indicate which of the following are accounted for: a reddening 
dependent correction (see Section l4.3l l. a per-chip zeropoint (dis- 
cussed in Section l2.3l l. a radial distortion term (Section l2.2l l. and 
a residual 2D spatial systematic correction (Section l3.2.3l l. 



Release Date E(B-V)' chip ZP radial spatial 



EDR 


Feb 2006 


no 


no 


no 


no 


DRl 


Jul 2006 


no 


no 


no 


no 


DR2 


Mar 2007 


yes 


no 


no 


no 


DR3 


Dec 2007 


yes 


yes 


no 


no 


DR4 


Jul 2008 


yes 


yes 


yes 


yes 



We demonstrate that for the vast majority of observa- 
tions with WFCAM, these tests are passed, and we quantify 
the conditions under which they are not. 

The paper is organised as follows: Section [2] summarises 
how the WFCAM data are calibrated within the pipeline; 
in Section |3] we examine the repeatability of WFCAM pho- 
tometry, quantify residual spatial systematics, and investi- 
gate the limits of the calibration in non-photometric condi- 
tions; in Section |4] we investigate the relationship between 
the 2MASS and WFCAM photometric systems and inves- 
tigate the effects of interstellar reddening; in Section [5] we 
investigate the extent of bias in the WFCAM calibration 
arising from the overestimation of flux in faint sources in 
2MASS; in Section [6] we compare the WFCAM system to a 
Vega system where AO stars should have zero colour in all 
passbands; in Section [7] we compare the WFCAM photom- 
etry with published data for faint standards measured with 
UKIRT; Finally, in Section [S] we derive the throughput of 
WFCAM. At the time of writing, UKIDSS has completed 
its fourth Data Release (DR4). 



2 METHOD FOR ROUTINE PHOTOMETRIC 
CALIBRATION OF WFCAM DATA 

All WFCAM data are processed in Cambridge with a 
pipeline developed by CASU and documented in Irwin et 
al. (2008). In this section we describe in detail the steps 
used to measure magnitudes and then apply a photometric 
calibration for each WFCAM image and subsequent cata- 
logue. After processing and calibration, the WFCAM data 
(in the form of FITS binary images and catalogues) are 
transferred to Edinburgh for ingestion into the WFCAM 
Science Archive (WSA, Hambly et al. 2008). The UKIDSS 
data have now seen a number of data releases from the WSA, 
and the photometric calibration has evolved to tackle the 
(increasingly smaller) corrections needed to meet the sur- 
vey goals; the corrections included in each release are sum- 
marised in Table [T] 

2.1 Key WFCAM elements 

WFCAM comprises four Rockwell Hawaii-II PACE arrays, 
with 2k X 2k pixels at 0.4 arcsec/pixel, giving a solid angle 
of 0.21 deg^ per exposure. A detailed description of WF- 
CAM can be found in Casali et al. (2007). A prerequisite for 
the calibration scheme is to ensure that all four detectors 



are calibrated to the same photometric system. The imple- 
mentation of the reduction pipeline needs to ensure that the 
different intrinsic detector quantum efficiencies (and pos- 
sibly their different colour dependence) are correctly taken 
into account. Thus much of the calibration is performed per- 
detector. 

There are five broadband filters designed for WFCAM, 
the ZYJHK filters, which are described in detail by Hewett 
et al. (2006) and characterized using synthetic photometry. 
The JHK filters possess transmission profiles specified to 
reproduce as closely as possible the widely adopted MKO 
system (Tokunaga, Simons & Vacca 2005). The Z and Y 
filters cover 0.84-0.93/^m and 0.97-1.07/im respectively. 

2.2 Measurement of instrumental magnitudes 

The source extraction software (Irwin et al. 2008) measures 
an array of background-subtracted aperture fiuxes for each 
detected source, using 13 soft-edged circular apertures of 
radius r/2, r/v'2, r, V^r, 2r ... up to 12r, where r — 1 
arcsecond. A soft-edged aperture divides the flux in pixels 
lying across the aperture boundary in proportion to the pixel 
area enclosed. In this paper we only consider photometry 
derived from fluxes measured within an aperture of radius 
1 arcsecond. However, all the apertures of selected isolated 
bright stars are used to determine the curve-of-growth of 
the aperture fluxes, i.e. the enclosed counts as a function of 
radius. This curve of growth is used to measure the point 
spread function (PSF) aperture correction for point sources 
for each detector, for each aperture (up to and including 
4r, which includes typically ~99per cent, or more, of the 
total stellar flux). Irwin et al. (2008) find that this method 
derives aperture corrections which contribute < 1 per cent 
to the overall photometry error budget 

A further correction should be applied to the source 
fiux to account for the non-negligible field distortion in WF- 
CAM, described in detail in Irwin et al. (2008). The astro- 
metric distortion is radial and leads to an increase in pixel 
area by around 1.2 per cent compared to the centre of the 
field of view. Standard image processing techniques assume 
a uniform pixel scale, and that a correctly reduced image will 
have a flat background. For WFCAM's variable pixel scale, 
this is actually incorrect and one would expect to see an 
increase in the sky counts per pixel at large off-axis angles, 
while the total number of counts detected from a star would 
be independent of its position on the array. The flatflelding 
of an image therefore introduces a systematic error into the 
photometry of sources towards the edge of the fleld of view. 
The corrected flux /cor, where / is defined as the aperture 
corrected count-rate in ADU per second of the source above 
background, is simply: 

/cor = ,f/(l + 3fc3r')(l-Kfc3r') (1) 

Where kg (with units of radian/radian^) is the coeffi- 
cient of the third order polynomial term in the radial distor- 
tion equation (Irwin et al. 2008) and is found to be slowly 
wavelength dependent with a value of -50.0 in the iif-band 
(called PV2_3 in the FITS headers). The instrumental mag- 
nitude is then 

= -2.5logio{fcoi) (2) 
This correction has not been applied for WSA releases 
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Figure 1. Photometric radial distortion (magnitudes) in WF- 
CAM as a function of offaxis angle (arcmin). There are curves for 
all five broadband filters, ZYJHK from bottom to top. The WF- 
CAM detectors are separated by 94 per cent of their active area 
in a square pattern (see Fig. |4]l , thus the vertical lines represent 
the location of the detector corners. 



DR1-DR3 (Hambly et al. 2008), but is included for DR4 and 
subsequent releases. Fig. [T] plots the radial distortion term 
(/cor// converted to magnitudes) as a function of off-axis 
angle for the WFCAM filters. 

2.3 Calibration of the photometry 

The data are firstly flatfielded using twilight flatfields (which 
are updated monthly), and an initial gain correction is ap- 
plied to place all four detectors on a common system, to first 
order. The per-detector magnitude ZP is then derived for 
each frame from measurements of stars in the 2MASS point 
source catalogue (PSC) that fall within the same frame. 
Thus the calibration stars are measured at the same time 
as the target sources, taking us away from more traditional 
calibration schemes, whereby standard star observations are 
interspersed with target observations. 

We assume that there exists a simple linear relation 
between the stellar 2MASS and WFCAM colours, e.g. 
Jw — J2 oc J2 — H2- In a Vega-based photometric sys- 
tem, this relation should pass through (0,0), i.e. for an AO 
star Z — Y = J — H = K, irrespective of the filter system in 
use. For each star in 2MASS observed with WFCAM, the 
pipeline derives a ZP (at airmass unity) from 

ZP = m2 + C( J2 - H2 or J2 - K2) -m, + k{x - 1) (3) 

where nii is the aperture corrected instrumental magni- 
tude (derived above). C are the colour coefficients for each 
passband and have been solved for by combining data from 
many nights (see Section |4]). fc is a default value for the ex- 
tinction (0.05 in all filters, see below) and x is the airmass. 



The 2MASS sources used for the calibration of WF- 
CAM are selected to have extinction-corrected colour 0.0 < 
J2 - K2 < 1.0 with a 2MASS signal-to-noise ratio > 10 
in each filter. If fewer than 25 2MASS sources are found 
within the field-of-view of the detector, then the colour cut 
is not applied. The WFCAM astrometric calibration is also 
derived from 2MASS and for both systems the astrometric 
rms accuracy per source is < lOOmas (Irwin et al. 2008) . The 
maximum allowed separation for a 2MASS- WFCAM source 
match is 1 arcsecond. For all WFCAM observations up until 
2008 March 28th, the median number of 2MASS stars falling 
within a WFCAM detector, which match our colour and 
signal-to-noise requirements, is ~ 200. Only around 0.1 per 
cent of WFCAM pointings have fewer than 25 2MASS stars 
covering a single detector. For the UKIDSS Large Area Sur- 
vey (LAS) alone, the median number of 2MASS standards 
per detector is ~ 100, and again, around 0.1 per cent have 
fewer than 25. 

For a single pointing, for each detector, the zeropoint 
is then derived as the median of all the per-star zeropoint 
values. A single photometric zeropoint for the pointing, 
MAGZPT, is calculated as the median of the detector zero- 
points over all 4 detectors. The associated error, MAGZRR, 
is 1.48 X the median absolute deviation of the detector ze- 
ropoints around the median. This error therefore comprises 
several components: the intrinsic errors in the 2MASS pho- 
tometry, the error in the conversion from 2MASS to the 
WFCAM filter, and then any residual systematic offsets be- 
tween the detectors (see Section [3.2.2p . In the Z-band, the 
MAGZRR errors are the largest, typically about twice that 
for the JHK fihers. 

It should be noted that we do not derive any atmo- 
spheric extinction terms on a given night with WFCAM. 
Rather, the value of MAGZPT derived above incorporates 
an instantaneous measure of extinction at the observed air- 
mass. The photometric calibration of a field therefore in- 
cludes no error from this assumption. However the derived 
zeropoint for airmass unity will include a small error (be- 
cause we assume the extinction is 0.05 magnitudes per air- 
mass for all filters). Leggett et al. (2006) find the extinction 
at Mauna Kea to be kj = 0.047, = 0.029, fcjf = 0.052 
with standard deviations between 0.2 and 0.3. For a typ- 
ical WFCAM frame, observed at an airmass~ 1.3, an ex- 
tinction which differs from our assumed value by 0.03 mag- 
nitude/airmass will lead to a 0.01 magnitude error in the 
value of MAGZPT. The value of MAGZPT over time can 
be used to investigate the long term sensitivity of WFCAM 
due to, for example, the accumulation of dust on the optical 
surfaces, and seasonal variations in extinction. 

2.4 Detector offsets 

A final stage to the photometric calibration takes account 
of systematic differences between the four detectors, mea- 
sured on a monthly basis. The object catalogues associated 
with each science product frame are used in the pipeline 
calibration process to compute a single per pointing overall 
zeropoint for the array of 4 detectors, using the colour equa- 
tions specified in Section |4l The residuals from all 2MASS 
stars used in the frame zeropoint determination (i.e. J, H, 
K signahnoise > 10:1) are also computed on a per point- 
ing basis together with their standard coordinate location 
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(C) v) with respect to the tangent point of the telescope op- 
tical axis (see e.g. Irwin et al. 2008). We use standard co- 
ordinates since these are independent of the degree of inter- 
leaving and dithering, and hence pixel scales, used in con- 
structing the science product images and catalogues. These 
residuals are then partitioned by month and stacked. The 
monthly timescale corresponds with the changover of the 
master flat field frames, with which we anticipate some of the 
corrections are correlated. From UKIDSS DR3, all WFCAM 
photometric calibration is determined and applied monthly. 
The details of this correction process are discussed in Sec- 
tion [3311] The result is to generate zeropoints per detector 
(the value of MAGZPT is updated for each detector) . 

The pipeline also estimates the nightly zeropoint 
(NIGHTZPT) and an associated error (NIGHTZRR) which 
can be used to gauge the photometricity of a night. 
NIGHTZPT is simply the median of all ZPs measured within 
the night, and NIGHTZRR is a measure of the scatter in 
NIGHTZPT. 



3 REPEATABILITY OF THE PHOTOMETRY 

3.1 Overlap analysis from the UKIDSS Large 
Area Survey 

We can use repeat observations of sources in UKIDSS to 
investigate the photometric accuracy and homogeneity of 
the survey and to compare with the pipeline errors. The 
pipeline errors are calculated from the source counts and lo- 
cal background only, and do not include corrections for bad 
pixels and other detector artifacts. In addition, they delib- 
erately do not include any contribution from the calibration 
procedure itself (i.e. the inherent uncertainties in the cal- 
ibrators) or additional corrections for residual systematic 
effects in the calibration. By comparing multiple observa- 
tions of sources near the edges of the detectors we are more 
sensitive to systematics from spatial effects, e.g. scattered 
light in the Hatfield, a variable point spread function, and so 
on (see Section I3.2|l . Hence this analysis demonstrates the 
worst case, and should help put upper limits on the error 
contribution from systematic contributions. 

Fig. [2] shows the rms diagrams for repeat measure- 
ments, for all point sources (classed as stellar in both ob- 
servations) in the DR3 release of the LAS. Repeat observa- 
tions arise from the small overlap of detector edges between 
the four pointings which go towards constructing a tile, as 
well as from observations of adjacent tiles. The rms for each 
magnitude bin is calculated from the standard deviation of 
the differences between the two measurements from all the 
repeats, aAmag/V^- The data used in Fig. [2]were observed 
in photometric conditions, with a frame photometric zero- 
point within 0.05 magnitudes of the median value for all 
observations made with the filter (ctzp = 0.03 magnitudes 
in all pass bands). We exclude sources at the very edges of 
the array, within 10 arcseconds of the detector edges. 

For a photometric calibration good to 2 per cent, then 
the median magnitude difference should be zero, and the 
rms of the distribution for the brighter stars should be 
(j/^ = 0.02. 

At the bright end, the rms (black crosses in Fig. [2]) 
is very close to the survey goal of 2 per cent. Unsurpris- 
ingly, the pipeline estimated random noise errors (red dots 
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Figure 2. rms diagram for repeat observations of sources in the 
DR3 dataset for the LAS (selected to be stellar in both obser- 
vations) for the YJHK filters (top to bottom), showing results 
for data taken in photometric conditions. Photometric conditions 
are defined as having a zeropoint within 0.05 magnitudes of the 
survey median value. Black crosses are 1.48 X the MAD (median 
absolute deviation) of the data in magnitude slices, i.e. equiva- 
lent to the Gaussian sigma of the magnitude bin. Red circles are 
the median pipeline errors for sources in the bin as reported in 
the WFCAM Science Archive for DR3. The blue horizontal lines 
at 0.02 magnitudes represent the target 2 per cent photometric 
reliability. 



in Fig. ^ significantly underpredict the observed distribu- 
tion at the bright end, indicating that systematic errors, 
dominate here at the 2 per cent level. Figure |3] combines 
measurements for all filters and compares the estimated ran- 
dom photon noise errors, E, derived by the pipeline to the 
measured errors M (the rms described above), which in- 
clude both systematic calibration errors and random com- 
ponents due to photon noise. The points have been fitted 
with a simple relation of the form = cE^ + where 
the systematic component s = 0.021 ± 0.001, the constant 
of proportionality c = 1.082 ± 0.014. In principle these can 
be used to update the default pipeline error estimates, how- 
ever, we suspect that c is slightly greater than unity due 
to a combination of factors relating to edge effects on the 
detectors and noise covariance arising from inter-pixel ca- 
pacitance within the detector (Irwin et al. 2008), which the 
pipeline error estimates do not allow for. 
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Figure 3. A plot of measured errors (Af ) against estimated errors 
(_E) derived from the rms distributions for LAS overlaps. The 
solid line is a linear fit to M-^ versus and is described in the 
text. 



3.2 Spatial systematics 

3.2.1 Measurement 



As described in Section 12.31 we store the standard coordi- 
nates rj) and WFCAM magnitude for each observation of 
a 2MASS star. For each month of data and for each filter we 
then compute a per star zeropoint relative to the previously 
computed overall zeropoint for each science product image 
for the entire array. These differences are then stacked within 
a fixed standard coordinate grid, of cell size ~ 1.2 x 1.2 ar- 
cmin covering the complete array of four detectors, making 
use of selected data from each month of observations. Only 
science products taken in photometric conditions are used 
in this process to minimise the effects of non-photometric 
residual structure. The average offsets in the grid of values 
are then filtered (smoothed) using a combination of a two- 
dimensional 3-pixel bimedian and bilinear filter. This latter 
step ensures a smooth variation of values over the grid by ef- 
fectively correlating the corrections on a scale of neighbour- 
ing grid points. This is necessary since, even with a month 
of data, the number of points per cell is still only typically 
25-100. The rms noise per cell from the 2MASS errors alone 
is ~lper cent. With smoothing, this noise drops to the few 
milli-mag level which is a negligible extra error with respect 
to the derived corrections. 

The residuals for each detector are then grouped and the 
median correction calculated. These median factors define 
the detector-level zeropoint corrections for each passband. 
The final stage is to apply the detector-level corrections and 
compute the residual systematics which are recorded as a 
correction table, and diagnostic plot. An example of the 
latter is shown in Fig. |4] and the possible origins of these 
patterns are discussed below. The correction tables are cur- 
rently available via the CASU web page^. 



^ http:/ /casu.ast. cam. ac.uk/surveys-projects/wfcam/data- 
processing/illumination-corrections 
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Figure 4. An example of the remaining monthly average spatial 
magnitude residuals for WFCAM J-band observations in 2005 
September after applying passband-dependent individual detec- 
tor zeropoint corrections. The residuals are shown in standard 
coordinates with respect to the tangent point of the observations, 
i.e. E is to to the left and N is to the top. 



3.2.2 Detector-to-detector offsets 

The detector zeropoint level corrections show some inter- 
esting effects when data from five WFCAM semesters are 
examined (see Fig. (Sjl . The small random jumps in the tem- 
poral sequences are most likely due to detector- level DC 
offsets caused by flatfield pedestal effects (these can also be 
seen in the sky corrections but in this case they are generally 
additive and do not impact object photometry). Addition- 
ally, « zblper cent, more constant trends are seen in the 
offsets. These are correlated within a passband but are dif- 
ferent between bands. 

We wondered if small QE variations between the detec- 
tors could give rise to slightly different colour equations, re- 
quired to relate the natural detector-|-filter-|-telescope pass- 
band to the 2MASS standards. The overall passband gain 
differences between the detectors are corrected at the flat- 
fielding stage but because the twilight sky used in the flats 
is generally a different colour from the majority of objects, 
residual colour equation differences between the detectors 
could be manifest as a small zeropoint offset. Analysis of 
WFCAM images shows that the twilight sky changes colour 
rather rapidly as the sky darkens, ranging from fairly neu- 
tral {J-K ~ 0.5), through extremely blue {J - K ~ -2.0), 
to the rather red dark sky {J — K ^ 2.5). Modelling of the 
twilight sky and typical calibration stars (using a range of 
synthetic spectra as input) suggest that the detectors would 
have to exhibit large differences in QE, at the 10 per cent 
level, across individual passbands to explain the observed 
effect. Hewett et al. (2006) found that the QE of a typi- 
cal Rockwell Hawaii-II detector increases approximately lin- 
early with wavelength, by 8 per cent over the entire WFCAM 
0.8-2.2 micron range. However we currently lack measure- 
ments of the actual WFCAM detector QEs. For the time 
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exclude observations made in the first 4 months of WF- 
CAM operations, when adjustments to the focal plane geom- 
etry were still being made, the average systematic difference 
between 2 and 4 arcsec diameter stellar fluxes is negligible 
(<iper cent) over the entire array. 
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Figure 5. An overview of the computed detector zeropoint cor- 
rections as a function of observation date for all of the WFCAM 
broadband filters: Z - black; Y - cyan; J - blue; H - green; K - 
red. The patterns broadly split into a constant offset component, 
and occasional random jumps in level believed to be caused by 
DC pedestal level offsets in the derived fiatfield frames. The nu- 
meric label on each panel refers to the detector with #1 SW, #2 
SE, #3 NE, #4 NW as seen on the sky 



being, the cause of the detector-to-detector photometric off- 
sets remains an open question. 

3.2.3 Residual spatial systematics and their possible 
causes 

After removing the detector- level differences, generic recur- 
ring patterns are visible in the final spatial correction, even 
across different passbands. This strongly suggests a common 
underlying cause and the three most likely contributors are: 

(i) spatially dependent PSF corrections - the residual 
maps are made from aperture-corrected 2 arcsec diameter 
flux measurements. 

(ii) low-level non-linearity in the detectors; 

(iii) scattered light in the camera which would negate the 
implicit assumption of, on average, uniform illumination of 
the field; 

We have ruled out spatial PSF distortions as a signif- 
icant contributor by comparing aperture-corrected 4 arcsec 
diameter fluxes with the 2 arcsec diameter measures. If we 



3.2.4 Residual non-linearities? 

The spatial magnitude residuals are generally larger toward 
the edges of the chips. To test the derived lookup tables, 
we repeated the LAS overlap analysis (Section l3.1|l . but this 
time after correcting the source photometry. To apply the 
corrections we used bilinear interpolation to compute the 
offset for each source in the standard coordinate system de- 
scribed above. The correction is applied additively. 

The resulting rms-diagram (Fig. [6]) shows a significant 
improvement compared to the version without spatial cor- 
rection (Fig. [5}. The median values of rms for bright stars 
(taking the mag=13.5 bin in all filters), in each passband 
are typically 0.002-0.004 magnitudes lower (Table [l]). 

We can try to gain some insight into the cause of the 
spatial systematics by looking at the dependence of the cor- 
rection on magnitude. This would show up any residual non- 
linearity effects for example. In Fig.[7]we plot the quadrature 
difference between the spatially corrected and non-corrected 
rms as a function of magnitude, i.e. the error contribution 
from the spatial systematics. We show the diagram for stars 
as well as galaxies, as any non-linearity effects would be less 
pronounced in spatially extended sources. Fig. [7] shows that 
the spatial systematics amount to a correction of about 1 per 
cent in all filters with no dependence on magnitude. We con- 
clude that the WFCAM system shows no significant residual 
non-linearity. Our conclusion is consistent with the analysis 
of Irwin et al. (2008) who show WFCAM is linear to at least 
1 per cent. 



3.2.5 Scattered light as a cause of the spatial systematics? 

Fig. [5] shows the correlation between normalised fiatfield 
level and the spatial magnitude residuals for J-band data 
analysed from 2005 September. The detectors with the 
largest and smallest fiatfield variations (#1 and #4) show 
no effects at the ^1 per cent level; whereas detectors #2 and 
#3 show a clear correlation between residual magnitude spa- 
tial systematics and fiatfield level. The spatial variation in 
the fiatfield intensity is dominated by changes in sensitiv- 
ity across the detectors (and vignetting within the optical 
train) and reaches a factor two for the worst detector (Ir- 
win et al. 2008). The spatial systematics seen in the stacked 
2MASS photometric residuals represent a 1 per cent error 
contribution, and it seems most likely to arise from (as yet 
unquantified) scattered light within the WFCAM system, 
probably present in all observations. 



3.3 Repeat observations of standard fields 

Standard star fields are observed every night with WFCAM, 
and provide another means of testing the repeatability of the 
photometric calibration. For this analysis we made use of 
data calibrated for UKIDSS DR3. Twenty-seven standard 
reference fields were selected and, for each field, reference 
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Figure 6. rms diagram for repeat observations of sources in the 
DR3 dataset for the LAS (selected to be stellar in both obser- 
vations) for the YJHK filters (top to bottom), after correcting 
for the spatial systematics as described in the text. Data are for 
photometric conditions only. Black crosses are 1.48 X the MAD 
(median absolute deviation) of the data in magnitude slices, i.e. 
equivalent to the Gaussian sigma of the magnitude bin. Red cir- 
cles are the median pipeline errors for sources in the bin as re- 
ported in the WSA for DR3. The blue horizontal lines at 0.02 
magnitudes represent the target 2 per cent photometric reliabil- 
ity. 
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Figure 7. Arms is equal to ~ '""^^spatiai) plotted 

against magnitude for the Y (cyan triangle), J (blue diamond), 
H (green filled circle) and K (red plus) filters. In the upper panel 
we plot the data for objects classified as stellar, while the lower 
panel is for extended sources, rms is the rms for photometry in 
the DR3 release of the LAS (as described in Fig. |2]l. rmSgpatiai 
is the rms for the same data but with photometry corrected for 
spatial systematics as described in the text. 
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stars were chosen to have a minimum of 100 detected counts 
in a 2 arcsecond diameter aperture. Fields typically contain 
between 500 and 10,000 suitable reference stars depending 
on Galactic latitude. The reference stars were then matched 
against subsequent observations of the same fields taken in 
good conditions (i.e. the frames pass the following criteria: 
seeing < 1.52 arcseconds, airmass < 1.8, NIGHTZRR < 
0.05, MAGZRR < 0.05, median source ellipticity < 0.2). 
On average, the reference fields have between 20 (in the Z- 
band) and 50 (in the A'-band) observations which match our 
criteria. The smaller number of selected Z-band observations 
is due to the larger error in the frame zeropoints ( Section 

For the repeat observations of the standard fields, the 
same 2MASS calibrators are used in each observation, thus 
we should expect to see an improvement in rms. The ac- 
tual median values for rmssTD (for the 13.5 magnitude bin) 
range from 0.011-0.015 magnitudes for bright stars, which 



0.6 



1.2 



Normalised flatfleld leve 



Figure 8. The correlation between normalised fiatfield level 
and the spatial magnitude residuals for J-band data from 2005 
September. Detectors #1 (black) and #4 (blue) show no strong 
correlation, however detectors #2 (red) and #3 show a well- 
defined correlation. 



are 0.004-0.006 magnitudes lower than the spatially cor- 
rected LAS rmsspatiai- 
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Table 2. rms at mag=13.5 for the WFCAM filters derived from 
repeat observations of sources in the LAS and standard star fields. 
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Figure 9. rms as a function of slty extinction for stars with 
magitudes in the range 12-15. The sample is drawn from the LAS 
DR3 and the rms values for each bin are computed from the 
magnitude differences between repeat observations of the same 
star, see Fig.[2]for details. At least one of the frames in each pair 
must have an extinction value in the range plotted (except for 
the best weather bins, which must have both frames taken with a 
ZP within 0.05 mags of nominal). The data are plotted by filter, 
YJHK from top to bottom. 



3.4 Calibration performance for data taken in 
non-photometric conditions 

The measured zeropoint for each frame depends directly on 
the transmission of the atmosphere. The presence of clouds 
leads to a reduction in the derived zeropoint. In Fig. |9]we 
investigate the rms for stars measured in different condi- 
tions to investigate whether there is a specific sky trans- 
parency beyond which the photometric conditions become 
unacceptable. We again compute the rms for magnitude dif- 
ferences between repeat measurements, but this time over a 
larger magnitude range, m=12-15 in all filters (it's clear 
from Fig. [5] that the rms in this magnitude range is largely 
dominated by systematic errors over photon counting). The 
data are plotted as a function of sky extinction as measured 
by the difference of the frame zeropoint from the median 
zeropoint for that filter for at least one of the frames (ex- 
cept for the 'clear' sky case, where both frames must have 
an extinction < O.OSmag. 

Fig. |9] demonstrates that, even in non-photometric con- 
ditions, where 10-20 per cent of the stellar flux is lost to 
cloud, the calibration is still very good, and still meets the 
2 per cent calibration requirements. 



3.5 Photometric repeatability: summary 

In Table[2]we tabulate the rms at mag=13.5, for each filter, 
for the different samples considered in the sections above. 
Examining the first three columns, the largest calibration 
improvement comes from the computation of a per-detector 









LAS 




STD 


Filter 


rmSDR2 




'•'""spatial 


'•mScloudy 


rmssTD 


Zw 










0.015 




0.023 


0.021 


0.019 


0.021 


0.013 


J„ 


0.021 


0.018 


0.016 


0.022 


0.012 




0.022 


0.017 


0.015 


0.021 


0.011 




0.023 


0.020 


0.016 


0.020 


0.011 



zeropoint, as included in the change to WFCAM processing 
for DR3 from DR2. The overlap analysis of the LAS DR3 
data shows that the calibration requirements for WFCAM 
are met, and that the photometry is repeatable to ~2per 
cent, even for sources which are close to the edge of the 
field of view, where residual systematics are larger. Investi- 
gation of these systematics suggests that they are most likely 
caused by scattered light present in the WFCAM flatfields, 
which contributes about 1 per cent to the overall photomet- 
ric error budget at large off-axis angles. Application of our 
corrections for these systematic residuals reduces the WF- 
CAM photometric error to only ~1.5per cent in J, H and 
K. 

Repeat measurements of standard fields, where the 
sources are uniformly distributed across the field of view, 
and the 2MASS calibrators do not change yield a photo- 
metric repeatability at the level of 1-1.5 per cent. 

Even observations taken through thin cloud (where 
transmission is > 80 per cent) can be calibrated to ~2per 
cent using 2MASS stellar sources as secondary standards. 



4 COMPARISON BETWEEN WFCAM AND 
2MASS PHOTOMETRIC SYSTEMS 

The colour equations determined for regions of low redden- 
ing are listed below. They were derived from fits to selected 
data taken in the first year of WFCAM observations (Irwin 
et al. 2008). The 0.03 magnitude offset for the ff-band is dis- 
cussed below. The F-band offset was not included prior to 
DR3, and is discussed in Section r6.2l In this Section we red- 
erive the 2MASS-WFCAM tranformations, and investigate 
their behaviour in regions of low and high Galactic latitude. 
We establish the reliability of the photometry for fields con- 
taining relatively small numbers of 2MASS calibrators. We 
also quantify the effects of Galactic extinction on the cali- 
bration, and derive additional corrections for Equations|3HHl 
listed in Equations l9l ll3l For DR2 onwards, it is this full set 
of equations that are applied to the 2MASS secondary stan- 
dards to enable derivation of the photometric zeropoint. 



= J2 +0.950(J2-//2) 

Yw = J2 + 0.500( J2 - Hi) + 0.080 

Jw = J2 - 0.065(J2 - H2) 

H^ = H2 + 0.070(J2 -H2)- 0.030 
= K2+0.010(J2~ K2) 



(4) 
(5) 
(6) 
(7) 
(8) 
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Table 3. Table of best fit colour terms for WFCAM-2MASS 
colour as a function of J2 — H2 (or J2 — K2 in the case of the K 
filter) as shown in Fig. 1101 split by region (see text for details). 
The first column shows the best fit colour terms for data taken 
in the second half of 2005 (semester 05B) and can be directly 
compared to the colour equations listed in Equations |4] to |8] 





Semester 




Region 




Colour 


05B 


A 


B 


C 


Zv, — J2 


1.02 


1.04 


0.92 


1.01 


yw - J2 


0.55 


0.51 


0.55 


0.56 


Jvi — J2 


-0.05 


-0.04 


-0.06 


-0.05 


Hv; — H2 


0.08 


0.04 


0.07 


0.07 


K^-K2 


0.02 


0.01 


0.01 


0.01 



4.1 The 2MASS colour-terms 

Fig. [To] shows data derived from observations taken in the 
second half of 2005. We plot the difference between the cal- 
ibrated WFCAM photometry and the 2MASS photometry 
for some 200,000 measurements across five passbands. We 
use these data to rederive the colour-terms to compare to 
those in Equations [IHll The data are selected to avoid re- 
gions of high interstellar reddening, E{B — V)0 

The data are fit with a linear relation (with 3cr clipping), 
including only stars with colours in the range 0.0 < J2 — 
K2 < 1.0. The slopes of the fits for are shown in Fig. llOl and 
can be seen to be in good agreement with Equations |4]to|8l 

For the JHK passbands, the linear conversions work 
surprisingly well, given the large spread in Figure 1101 For 
the //-band, an additional 3 per cent offset is required to 
bring the WFCAM photometry onto the same system as the 
UKIRT faint standards (as originally published by Hawar- 
den et al. 2001). The //-band offset is discussed in more 
detail in Section [7] 

For the Z and Y passbands, significant (and indeed de- 
generate) residuals are seen for stars with J2 — H2 ~ 0.5 — 
0.7. This is most easily understood by appealing to the two- 
colour J — H vs. H — K diagram for dwarf stars, e.g. Fig. 2 
in Cruz et al. (2003). Stars with spectral types later than 
MO become progressively bluer in J ~H, but redder in Z ~ J 
and Y — J. Thus the same colour transformation is applied 
to stars with the same J — H colours but very different 
temperatures and therefore Z — J colours, resulting in the 
structures seen in Fig. [10] 

The formal errors on all the fits are extremely small due 
to the large numbers of sources present in each sample, and 
are not particularly helpful when considering that they do 
not account for systematic effects (e.g. population). Below 
we consider two issues; (1) how robust is the determina- 
tion of zeropoint assuming the colour equations are correct 
for randomly selected small samples of stars, matched to 
the example of the Large Area Survey?; (2) are the colour 
equations significantly different in regions of high and low 
Galactic latitude? 



^ Throughout the paper we correct the E{B—V) from Schlegel et 
al. (1998) according to the prescription given by Bonifacio, Monai 
& Beers 2000, hence the use of the 'prime' on the formula. 




-0.2 |- 



0.0 0.2 0.4 0.6 0.8 

(J-H)2MASS 




-0.2 r 



0.0 0.2 0.4 0.6 0.8 1.0 

(J-K)2MASS 

Figure 10. A Hess diagram showing differences between WF- 
CAM and 2MASS photometry as a function of 2MASS ( J2 - H2 
or J2 — K2) colour for data taken in the second half of 2005. 
The data have been selected such that E(B — V)' < 0.2 and 
J2 — K2 < 2.0 to avoid the most significant effects of reddening. 
The best fit linear regressions to the unbinned data are overplot- 
tcd. 



4.2 Robustness of zeropoint determination 

For the determination of the zeropoint of a typical high- 
latitude field, the measured offset between the WFCAM in- 
strumental magnitudes and the 2MASS photometry will not 
be affected by these colour-dependent residuals, so long as 
a robust (i.e. median) offset is measured. Clearly the av- 
erage offset in the top two panels in Fig. [TD] would be af- 
fected by the non-linearity in the colour-transformations. If 
a WFCAM field were dominated by extremely red stars, or 
an unusual population of objects (for example a large and 
rich Globular Cluster), then a further colour-dependent ze- 
ropoint correction may be necessary. We have yet to identify 
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Table 4. Average Zeropoint offsets and their standard deviations 
for a comparison between the 2MASS and WFCAM photometry 
for 10000 trials with randomly chosen samples of 50 stars 



Filter 


Mean ZP Offset 


'^'offset 




-0.006 


0.017 




-0.001 


0.011 




0.003 


0.007 


Hw 


0.002 


0.008 




0.003 


0.010 



any fields that suffer from this effect. For example, examin- 
ing zero-point variation as a function of Galactic latitude 
and longitude reveals no compelling evidence for any vari- 
ation significanlty larger than 1 per cent. We anticipate re- 
turning to this issue in the future as more data are collected 
over larger areas of sky. 

The robustness of the calibration of WFCAM is most 
simply tested using a Monte-Carlo approach. We took a 
'worst-case' scenario, comprising a high-latitude field which 
contains relatively few 2MASS stars. For each trial we se- 
lected 50 stars at random (from the dataset used above), and 
derived the median offset between the 2MASS magnitudes 
corrected to the WFCAM system (using Equations [JHH]) and 
the WFCAM calibrated photometry. We performed 10,000 
trials, and summarise the results in Table U This method of 
calibration is found to be robust for all filters, with a stan- 
dard deviation of 1 per cent or lower in all filters, except the 
Z-band, for which the error is closer to 2 per cent. 



4-2.1 Colour equations at low and high Galactic latitude 

We defined three regions to compare: A (0 < / < 100, 
|6| < 10), B (0 < Z < 100, l&l > 20), C (180 < I < 250, 
|6| > 15), where I and b are Galactic longitude and latitude 
respectively. Region A thus contains a fairly large sample 
close to the Galactic plane. As with the sample examined 
above, we fit for the colour transformations between the 
2MASS and WFCAM filters. The results of the fitting are 
listed in Table [3] There appear to be some small differences 
between the regions, particularly in the Z-passband for re- 
gion B. However, bootstrap error estimates find an error in 
the slope of ±0.02 for regions B and C in this filter, so we 
do not believe these differences to be significant (above 2a). 
Note that the colour term for a pointing would need to be 
adjusted by 0.04 to change the zeropoint by 2 per cent for 
stars of median J2 — H2 ~ 0.5. 



4.3 The Effects of Galactic extinction 

We have attempted to mitigate against the effects of inter- 
stellar reddening on the WFCAM calibration by employing a 
colour restriction to the stars used in the calibration (applied 
after extinction correction) . At regions of low Galactic lati- 
tude, where reddening is high, the spectra of distant stars are 
strongly affected by intervening dust. The population mix 
of the WFCAM calibrators will change as intrinsically blue 
stars are moved into the applied colour selection, and red- 
der stars are moved beyond the colour bounds. Giant stars 
become significantly more common. Differences between the 
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Figure 11. Frame zeropoint minus nightly zeropoint (uncor- 
rected for reddening) versus E(B — V)' for the WFCAM broad- 
band filters. Only data from photometric nights are included. 
Data points which lie more than 0.1 magnitudes away from the fit 
are clipped. The blue line is the best fit relation to the extinction 
relation, while the red line is constant AZP = 0.0. 



WFCAM and 2MASS filter profiles become important, and 
the colour equations used to transform the 2MASS magni- 
tudes onto the WFCAM system could begin to break down. 
For the JHK filters, for a mixed stellar population, these 
differences are small, and we expect the effect on the calibra- 
tion to be similarly small. For the Z and Y filters, where we 
are extrapolating from the J2 magnitudes using the J2 — K2 
colour then we may expect to see significant reddening de- 
pendent offsets in the calibration. 

The following analysis makes use of data calibrated ac- 
cording to the prescription for DRl, which made no attempt 
to correct for reddening. In Fig. [TT] we plot the difference 
between the nightly averaged zeropoint and the zeropoints 
measured for individual fields taken on that same night, ver- 
sus E[B — V)' (corrected according to Bonifacio, Monai & 
Beers 2000) for the ZY JHK filters using data taken on 
photometric nights between 2005 April and 2007 May. The 
E(B — Vy value for each frame is computed as the average 
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Figure 12. The standard deviation of AZP (computed from the 
MAD) around the best fit relations described in Equations |5] to 
[T3las a function of E{B - V)' . The data are binned in E{B - V)' 
as indicated by the horizontal bars. 



-E(i3 — 1^)' of all stars contributing to the calibration of the 
detector. Data points which lie > 0.1 magnitudes outside 
the best fit slopes are clipped. 

A significant, E{B — V)' dependent, correction is clearly 
required for the Z and Y filters. The correction required 
for the J, H and K passbands is more than an order of 
magnitude smaller, for example ~3 per cent in J for a change 
in AE{B~Vy = 2.0. Fits to these data are shown in Fig.fTTI 
and listed in Equations refeqzpz- H3l The colour equations 
employed for UKIDSS DR2 onwards have been modified to 
include these additional E{B — V)' terms. 

ZP'{Z) = ZP{Z) - 0.390E{B - V)' (9) 

ZP'(Y) = ZP{Y) - 0.160E{B - V)' (10) 

ZP'{J) ^ ZP{J) -0.015E{B -Vy (11) 

ZP'{H) ^ ZP{H) -0m5E{B -Vy (12) 

ZP'{K) = ZP{K) -0m5E{B -Vy (13) 

There exists a paucity of measurements at high redden- 
ing for the Z and Y passbands, because they are not in- 
cluded in the Galactic Plane Survey. Nevertheless, we inves- 
tigate the scatter in the measured zeropoints as a function of 
E{B — V)' and try to estimate at what point the reddening 
becomes sufficiently large that the photometric calibration 
is compromised. This is attempted in Fig. [12] where we plot 
the standard deviation (computed from the median absolute 
deviation) of the data around the fits (Equations [9l to [T3)) as 
a function of E{B — Vy . azp is somewhat larger for the Z 
and Y bands at all values of E{B — V)' , as expected. For the 
JHK bands, the scatter in the measured ZP is seen to be 
very slowly rising as a function of E{B — Vy . Our conclusion 
is that the 2MASS-based calibration reaches the WFCAM 
goal for E{B - Vy < 2.0. For the Z-band, the WFCAM 
calibration goals are not met for E{B — Vy > 0.2, while for 



the F-band, the WFCAM calibration appears to be robust 
up to E{B - vy = 1.5. 

In summary, and more qualitatively, the zeropoints, and 
therefore the calibration, derived for highly reddened fields 
appear to be robust in the J-, H- and -ff- bands, but should 
be treated with some caution in the Z- and Y-bands. 



5 THE EFFECTS OF THE OVERESTIMATION 
OF FLUX FOR FAINT SOURCES IN 2MASS 

By calibrating WFCAM from a shallower survey, we open 
ourselves up to the possibility of a photometric bias affect- 
ing the results. All measured fiuxes in a survey are subject 
to uncertainties, and close to the detection threshold, where 
the uncertainties are large, this effect will result in a bias, 
whereby sources with negative fluctuations will not be de- 
tected, while those with positive fluctuations will be pushed 
to brighter magnitudes. Thus, on average, sources near to 
the detection threshold of a survey will have overestimated 
fluxes. This is generally known as the Eddington Bias (Ed- 
dington 1940). 

In order to investigate the impact of Eddington Bias on 
the WFCAM calibration, we make use of the 2MASS Cal- 
ibration Merged Point Source Information Tablei|f|. These 
tables contain the results of merging catalogues from the 
individual 2MASS calibration scans to give average mea- 
sured photometry and astrometry for all sources observed 
multiple times. We found 7 of the 2MASS calibration scans 
overlapped with WFCAM standard flelds with at least 5 
observations (though most had > 100). 

In Fig. [T3] we compare magnitudes from the 2MASS 
merged tables with observations of the same sources in WF- 
CAM. Note that the WFCAM magnitudes combine mea- 
surements of at least five observations for all objects. The 
WFCAM calibration is derived from 2MASS point sources 
with a signal-to-noise ratio > 10. As can be seen in Fig. 1131 
the Eddington Bias at these levels is < 2 per cent in all fil- 
ters, indicating that the error introduced into the WFCAM 
calibration is negligible. We also note for the J-band, there 
is a possible trend between magnitude and the magnitude 
difference, suggestive of a small non-linearity. The effect is 
less pronounced at H and probably non-existent at K. 



6 THE WFCAM ZYJHK SYSTEM 

WFCAM photometry is on a Vega system (see Hewett et 
al. 2006). The photometric zeropoints for aU the WFCAM 
fllters for each observed frame are derived by measuring the 
offsets between the 2MASS calibrators (now converted to 
the WFCAM system) and the observed stars, as described 
in Section [2] using the colour equations from Section [4] For 
the Z and Y passbands in particular, there are no ready 
calibrators available, and there is an underlying assumption 
that the 2MASS colours can be linearly extrapolated into the 
WFCAM system. In this section we investigate the WFCAM 



^ http:/ /www. ipac.caltech.edu/2mass/reIeases/allsky/doc/ 
seca6_l.html 
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Figure 13. Differences between measured WFCAM photometry for calibration fields and the same sources in the 2MASS Calibration 
Merged Point Source Information Tables (www.ipac.caltech.edu/2mass/releases/allsky/doc/seca6_l.html) plotted as a function of 2MASS 
magnitude for all filters, illustrating the effects of Eddington Bias. The lOo" magnitude thresholds for calibration star selection arc shown 
as a vertical line. 



photometric system in more detail using colour-colour dia- 
grams, and in particular concentrating on photometrically 
selected AO stars. 



6.1 Sources in the SDSS 

We define a sample of stellar sources detected in both the 
UKIDSS Large Area Survey (LAS) and SDSS. In addition 
we highlight stars with SDSS photometry consistent with 
a spectral type AO, i.e. having Vega colours close to zero 
{u — g,g — r,r — i,i — z all in the range -0.1 to O.lQ). This sam- 
ple is used to investigate the WFCAM photometry and check 
that the infrared colours are also close to zero. The J — H vs 
H — K diagram is shown in Fig. [T?] The median WFCAM 
colours for the candidate AO stars are J — H = 0.014 ±0.024 
and H — K — —0.008 ± 0.037 and are consistent with zero 
(with errors derived as the robust standard deviation of the 
data, a = 1.48 x MAD). In conclusion the stellar sequence is 
not offset in the JHK diagram. Synthetic WFCAM photom- 
etry is presented by Hewett et al. (2006) derived from the 
Bruzual-Persson-Gunn-Stryker (BPGS) spectrophotometric 
atlas. In Figure[l4]we plot the synthetic stellar sequence and 
we can see that there is a significant offset from the data. 
At this point it's not clear whether the offset should be at- 
tributed to either (a) the initial calibration of the BPGS 
spectrophotometry, (b) the synthetic colours generated by 
Hewett et al. (2006), or (c) systematic shifts present in the 
WFCAM calibration. 

We can explore this in a little more detail for each 
of the LAS filters in turn. We use the u — K colour to 
extend the baseline on the x-axis, and plot all WFCAM 
colour combinations in Fig. [TS] for the YJHK filters. Off- 
sets are measured from straight line fits to the blue stars. 

^ We convert all SDSS photometry onto a VEGA system, apply- 
ing the offsets from AB magnitudes listed in Hewett et al. 2006 
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Figure 14. The J — H vs H — K diagram for stellar sources in the 
UKIDSS LAS with counterparts in SDSS. Objects are selected to 
be stellar in all WFCAM filters and an overall stellar SDSS clas- 
sification. Measured photometric errors are < 0.05 magnitudes in 
each WFCAM passband. Sources with SDSS colours close to zero 
{u — g, g — r, r — i,i — z all in the range -0.1 to 0.1) are plotted 
with errors. Red crosses are synthetic photometry from Hewett 
et al. (2006) based on the BPGS atlas (see text for details). 



The slopes of the fits are very shallow, and thus any sys- 
tematic error in the u'- or K-hand photometry will not 
significantly affect the measured offset. At u ~ K = 0.0, 
we find Ayj = -0.075 ± 0.007, Ayh = -0.075 ± 0.018, 
Ayk = -0.117 ± 0.018, AjH = -0.009 ± 0.012, Ajh = 
-0.009 ± 0.012, Ahk = -0.022 ± 0.024, indicating an ap- 
preciable offset for the K-band, and little or no offset in the 
JHK fihers. 
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Figure 15. WFCAM colours for stellar sources in common with 
SDSS. Large symbols are selected to be blue in SDSS (u — g,g — 
r, r—i, i—z all in the range -0.1 to 0.1). Error bars are the standard 
deviations on binned median values. The best fit straight lines are 
also drawn. 



6.2 The y-band 

At the start of WFCAM observations, there was no reUable 
cahbration for Y , and as a consequence an estimate of the 
zeropoint was made for the early data releases. To examine 
the y-band calibration in more detail, we plot the y — J vs 
J — K diagram in Fig. 1161 The blue star sequence does not 
pass through (0,0), but is significantly shifted to the blue 
in y — J. We attribute the displacement to the offset in the 
y-band calibration in DR2. A maximum likelihood straight- 
line fit to the blue stars (including photometric errors, see 
Appendix) yields the intercept in y — J, Ay = —0.074 ± 
0.005. 

Relaxing the source selection such that there is no 
requirement for a counterpart in SDSS yields a slightly 
larger sample. From the LAS dataset alone we measure 
Ay = -0.080 ± 0.005. 

Fig. 1 16 1 shows that there is a non-linearity in the trans- 




J-K 



Figure 16. The DR2 Y — J ws J — K diagram for stellar sources 
common to the UKIDSS LAS and SDSS surveys. Objects are 
selected to be stellar in all WFCAM filters with a measured pho- 
tometric error < 0.05 magnitudes in each passband. SDSS sources 
are selected to have a stellar classifcation and a photometric er- 
ror < 0.1 in all filters. Sources in regions of significant reddening, 
E{B — V)' > 0.2 are excluded from the plot and from analysis. 
So-called 'blue stars' are plotted with a larger symbol. The best- 
fit straight line is shown in blue. An additional linear fit to the 
remaining stars is shown in red extending down to J — = 0.8. 



formation between the J — K and Y — J colours for blue 
stars. Hence the y-band offset is a natural consequence of 
adopting a linear extrapolation based on redder stars. 



6.3 The Z-band 

The only UKIDSS survey with significant Z-band coverage 
is the Galactic Clusters Survey (GCS), with essentially no 
overlap with SDSS at present. In Fig. [T7]we plot Y — J and 
Z — J against J — K, for stellar sources detected in the GCS 
to: (i) confirm the y-band offset discussed above, (ii) check 
the validity of using the GCS dataset to measure the offset, 
and (iii) investigate a possible offset in the Z-haxiA. 

The sample of blue stars is rather smaller than derived 
from the LAS sample, but the offset for the y-band is re- 
produced with Ayj = -0.095 ± 0.015. For the Z-band we 
measure an intercept of Azj = —0.008 ± 0.015, i.e. consis- 
tent with zero. 



6.4 Summary of offsets 

In summary, an offset to the WFCAM y-band zeropoint 
needs to be applied to bring it onto the Vega photometric 
system. We take the weighted mean of the various values 
derived above with respect to the WFCAM J-band and the 
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SDSS i-band, and find Ay = -0.078 ± 0.010. Consequently 
an offset of 0.080 has been applied for UKIDSS DR3 and 
subsequently as shown in Equation [5] 

For the ZJHK-hanAs, the data are consistent with no 
offset at the ±2 per cent level. 



Figure 18. Median Amag (WFCAM minus UFTI photometry) 
plotted as a function of J — K (UFTI) for the J, H and K filters 
(top to bottom). The dashed lines show the median measured 
offsets between the WFCAM and UFTI photometry, discussed in 
the text. 



7 COMPARISON WITH PUBLISHED JHK 
PHOTOMETRY IN THE MAUNA KEA 
OBSERVATORIES PHOTOMETRIC SYSTEM 

Every night that WFCAM is observing, a number of stan- 
dard fields are measured, typically every two hours (in the 
first year of operations this was an hourly procedure). The 
majority of the standard fields include stars selected from 
the list of UKIRT faint standards published in Hawarden 
et al. (200lfl These fields were initially chosen to act 
as the primary source of photometric calibration for WF- 
CAM, however, subsequent analysis (presented in this pa- 
per) has demonstrated that the 2MASS-based calibration 
offers significant advantages. Observations of standard stars 
currently continue to be performed at the telescope, to en- 
sure that a future calibration can be derived independently 
from 2MASS should the need arise. The standard observa- 
tions also enable monitoring of the WFCAM system by the 
reliable strategy of repeatedly looking at the same stars. 
Recently, Leggett et al. (2006) have revisited the 



^ http:/ /www. jach.hawaii.edu/UKIRT/astronomy/calib/phot_cal/ 



UKIRT faint standards, and present new measurements us- 
ing the UKIRT Fast Track Imager (UFTI). They find small 
offsets (at the few per cent level) with respect to the older 
data, as well as some evidence of magnitude-dependent non- 
linearity, and recommend adopting the newer photometry. 

In Fig. [18] we compare the DR3 WFCAM calibration 
with the recalibrated UFTI photometry of faint standard 
stars measured in the Mauna Kea Observatories system 
(Leggett et al. 2006). 

WFCAM measurements were selected from observa- 
tions taken on photometric nights during good conditions. 
There are typically 50-100 measurements per star in each 
filter. The data values plotted are the median mw — mupTi 
for the JHK filters, and error bars are the standard de- 
viation calculated as 1.48 x the median absolute deviation 
(MAD) of the WFCAM measurements. Only stars passing 
the following criteria are plotted and analysed: (i) they show 
no variability {rmsjHK <~ 0.03); (ii) they are fainter than 
JHK=11 to avoid possible saturation/non-linearity effects. 
Thus the sample is rather small, and does not include any 
objects with extreme colours. 

Most stars are not variable at the 1-2 per cent level, 
but we note a few exceptions: FS117 (b216-b9), FS143 (Ser- 
EC86), GSPC P259-C, which all show evidence for variabil- 
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ity at > 5 per cent within the timeframe of the observations 
(2005 April-2007 May). As in Leggett et al. (2006), we find 
FS144 (Ser-EC84) may be variable at the 3 per cent level. 
However, we find no evidence to suggest that FS116 (b216- 
b7) is variable at > f per cent. 

Over the limited colour range considered in Fig. 1181 we 
do not see evidence for a colour term between the Jw and 
JuFTi photometry. However, a larger sample of stars with 
both UFTI and WFCAM photometry, spanning a broader 
colour range, is required before a definitive statement can 
be made. The median measured offsets between the WF- 
CAM and UFTI photometry are: Sj = 0.005, Sh = 0.015, 
Sk = 0.017. WFCAM-UFTI is positive in the three filters, 
i.e. the WFCAM measurements are slightly fainter. Leggett 
et al. (2006) also measured small offsets between the WF- 
CAM and UFTI systems; and find a different offset in the 
H-hand, with — -ffuFTi ~ —0.015. We note that the 
Leggett et al. (2006) measurements are based on an ear- 
lier WFCAM calibration, and use a sample of stars without 
such a restrictive magnitude cut and comprising many fewer 
repeat observations. 

During the WFCAM Science Verification phase (2005 
April), a decision was taken to anchor the WFCAM calibra- 
tion to the Hawarden et al. (2001) UKIRT photometry. An 
offset of 0.03 magnitudes was found between the 2MASS- 
based WFCAM cahbration and the UKIRT system, thus 
requiring a correction to the colour equation (Equation [T]) 
for the H-hand. When making this decision, we were also 
aware that Cohen et al. (2003) argued for small offsets to 
be applied to the 2MASS H2 and K2 bands to bring them 
into agreement with their space-based infrared spectropho- 
tometric calibration system (the offset required for J2 is 
negligible). The suggested offsets, to be applied to the pub- 
lished 2MASS photometry, were -0.019 ± 0.007 {H2), and 
-h0.017 ± 0.005 {K2). 

At present, in the absence of a larger sample of observa- 
tions and more-accurate photometry, we conclude that the 
WFCAM photometric system is a Vega system (to within 
2 per cent) and is in agreement with all three of the Hawar- 
den et al. (2001), Cohen et al. (2003) and Leggett et al. 
(2006) systems to within 2 per cent. For the time being we 
choose to retain the 3 per cent offset in the conversion be- 
tween the 2MASS and WFCAM Jf-band photometry. 



8 EMPIRICAL DETERMINATION OF THE 
THROUGHPUT OF WFCAM 

8.1 Gain 

Measurements of the readout noise and gain from the 2005 
April SV data are given in Table [S] The readout noise was 
estimated from the difference between successive single ex- 
posure CDS 5 s dark frames after running the decurtaining 
algorithm to remove systematic artifacts. Note that the av- 
erage dark current is generally negligible and dominated by 
reset anomaly variations. The gain is the average of gains 
measured at three background levels 23k, 14k and 5k; the 
overall variation in measured gain was at the 0.05 level with 
no clear trends with background. 

Whilst measuring the gain from the dome fiat se- 
quences, an estimate of the inter-pixel capacitance via a ro- 
bust measure of the noise covariance matrix was also made. 



Table 5. Measured dark current, uncorrected gain and noise (in 
electrons and ADU) for the four WFCAM detectors. The gain is 
an overestimate, as described in the text. The corrected average 
gain for WFCAM is 4.31. 



Det 


Dark curr. 


Noise (ADU) 


Uncorr. gain 


Noise (e-) 


#1 


-2.4 


3.8 


4.84 


18.4 


#2 


3.7 


4.2 


4.87 


20.5 


#3 


-9.0 


4.2 


5.80 


24.4 


#4 


49.1 


4.4 


5.17 


22.7 



A sum of the noise covariance matrix close to 0,0 gives a 
value of 1.20 which implies that the total reduction in di- 
rectly measured noise variance (i.e. measured using a con- 
ventional method) is therefore also 1.20. We therefore pre- 
dict a 20 per cent overestimate of the gain, and hence a 
20 per cent overestimate of the QE coefficients. Such an over- 
estimate is as expected for the Rockwell Hawaii-II devices. 
The average gain of the WFCAM detectors is therefore 4.31. 



8.2 Throughput 

The total throughput of the system is relatively easy to mea- 
sure (but it is much harder to quantify where in the system 
the actual losses are occurring) . For the calculation we have 
assumed the following: The effective area of the UKIRT pri- 
mary mirror is 10. Sm'^ (i.e. outer diameter 3.802m with an 
inner diameter of 1.028m). No attempt has been made to al- 
low for the shadowing of the primary by the secondary, nor 
for any obscuration caused by the forward mounted camera 
itself. The spectrum for Vega is the same used by Hewett 
et al. (2006), originally from Bohlin & Gilhland (2004). The 
passband transmissions are from Hewett et al. (2006) but 
renormalized to give the relative throughput as a function 
of wavelength, rather than absolute transmission (i.e. the 
peak value is 1.0). 

Thus Table [7] gives the estimated number of photons 
that would be incident on the primary mirror, assuming 
no atmosphere, multiplied by the relative transmission of 
the filter in each band. The values in table [7] allow the 
calculation of the zeropoint of the system in each filter 
assuming no losses due to the atmosphere, telescope and 
the instrument. Comparison with the median measured ze- 
ropoints for DRl then give us the throughput for WF- 
CAM+UKIRT-l-atmosphere. 



9 CONCLUSIONS 

(i) The CASU pipeline photometric calibration of WF- 
CAM data, using 2MASS sources within each WFCAM 
pointing, has achieved a photometric accuracy better than 
2 per cent for the UKIDSS Large Area Survey released in 
DR3 (Data Release 3) and later, based on repeat measure- 
ments. 

(ii) By stacking spatially binned WFCAM data on 
monthly timescales we have shown that photometric cali- 
bration residuals are present at the 1 per cent level. We have 
derived and tested a method for removing these residuals, 
and show that the photometric calibration for the WFCAM 
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Table 6. Estimate of WFCAM throughput 



Filt N„(0'") Nc(0"^) ZPioo ZP (M) TPm 





3.7cl0 


8.5c9 


24.95 


22.77 


15% 


Vw 


S.lelO 


7.3e9 


24.77 


22.78 


18% 


t/w 


2.9el0 


6.8c9 


24.70 


22.97 


23% 


H„ 


2.8cl0 


6.6c9 


24.66 


23.22 


29% 




1.6el0 


3.6e9 


24.02 


22.55 


29% 



Table 7. Where Ni/ is the number of photons reaching the de- 
teetor in each passband for a zeroth magnitude star, Nc is the 
equivalent number of counts in the detector and ZPioo is the 
corresponding predicted zeropoint for each filter if there were no 
losses in the system. ZPti,] is the actual measured median zero- 
point for WFCAM in each filter, and TPm is therefore the ad- 
justed throughput of the system compared to a perfect instrument 
and atmosphere. 

JHK filters can reach an accuracy of 1.5 per cent. We at- 
tribute these residuals to scattered light in the WFCAM twi- 
light flatfield frames. These corrections have been applied to 
the UKIDSS DR4. 

(iii) Even observations taken through thin cloud (where 
transmission is > 80 per cent) can be calibrated to ~ 2 per 
cent using 2MASS calibrators. 

(iv) Monte- Carlo sampling the 2MASS stars observed by 
WFCAM, suggests that as long as ~ 50 calibrators fall 
within the fiold-of-vicw, then the WFCAM calibration is ro- 
bust to 1 per cent in Y JHK and 2 per cent in Z. 99 per cent 
of the sky has more than 50 2MASS calibrators in a single 
WFCAM pointing. 

(v) The WFCAM calibration incorporates colour restric- 
tions on secondary standards and modified colour equations 
to mitigate against the effects of Galactic extinction. To- 
wards regions of significant reddening, the calibration meets 
the 2 per cent requirement for E(B — V)' < 2.0 for the 
JHK bands, but should be treated with caution above 
E{B — Vy = 0.2 for the Z band and to a lesser extent 
for the Y band. 

(vi) We find that the calibration of WFCAM from the 
shallower 2MASS survey is free from the effects of Edding- 
ton Bias, thanks to a careful choice of a signal-to-noise cut 
(SNR= 10) for the 2MASS sources. 

(vii) We identified a small offset between the WFCAM 
Y photometry and an idealized Vega photometric system 
in the second UKIDSS Data Release (DR2). A correction 
of Ay = —0.08 magnitudes has now been included in the 
WFCAM calibration and applied from UKIDSS DR3, in the 
sense that the Y band sources are now 0.08 magnitues fainter 
than in previous releases. 

(viii) We find that the WFCAM photometric system is 
within 0.02 magnitudes of the UKIRT MKO system (Leggett 
et al. 2006) for the JHK filters, and is Vegarlike to within 
2 per cent (i.e. stars selected to have AO colours in SDSS 
have WFCAM colours <= 0.02 across all passbands). 
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APPENDIX A: STRUCTURAL ANALYSIS 

Wc use the blue end of the distribution of points culled from 
the GCS to illustrate the problem of finding an accurate 
estimate for the F-band offset. In this case, since the WF- 
CAM JjiiT-bands are already on the Vega system, this is 
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Figure Al. Locating tlie offset in the Y-band calibration for 
WFCAM using standard least-squares regression (green and blue 
lines) compared to the correct Maximum likcUhood technique (red 
line). Red filled points are excluded from all the fits. 



equivalent to finding the intercept on the Y — J axis when 
the J — K colour is zero. However, with significant errors on 
both axes, and possibly also an intrinsic spread in the stellar 
locus, conventional least-squares curve fitting methods can 
give seriously biased results. This is illustrated in Fig. 
where the obvious outliers (red points) have been excluded 
from the fits. In this example we have (arbitrarily) defined 
the X-axis to be the J — K measures and the y-axis to be 
the Y — J measures. The green dashed line shows the result 
of a standard variance-weighted least-squares regression of 
y on X, while the blue dashed line shows the same for a re- 
gression of X on y. The values of the supposedly equivalent 
slopes and intercepts are significantly different and the clear 
bias in both fits is simply the result of ignoring the duality of 
the error distribution. We note that an often-used shortcut 
for this type of situation is simply to "average" the results 
from the two regressions and use that for the solution. As 
it happens, if the magnitude of the errors on both axes are 
similar this often gives a result close to the optimum value. 

The correct procedure in this type of case is to explic- 
itly recognise the more complex error distribution involved 
by rephrasing the modelling problem in terms of the un- 
observed "true" parameters using the method of maximum 
likelihood. This technique is also known as structural anal- 
ysis {e.g. Kendall & Stuart 1979). 

In general terms the underlying model is of the form 
Yi = f{Xi I 9), whereas what we observe are not these values 
but rather the quantities Xi = Xi -\- Sxi and yi = Yi + Syi 
where Sxt and Syi are independent errors with variance ct^. 
, CTy. which may include intrinsic spread in the model as 
well as observational and calibration errors. If these errors 
are Gaussian distributed then the likelihood L of observing 
the independent pairs of values Xi,yi is given by 



L^Yl P{x,,y, j 9) 



(Al) 



where 9 represent the model parameters; and the log- 
likelihood is therefore given by 




-0.15 -0.1 -0.05 
Intercept 



Figure A2. Log-likelihood contours starting at l-cr and spaced 
by intervals of l-cr for the maximum likelihood model fit shown 
for the data presented in figure El 



For unknown errors the problem is insoluble. For known 
errors 

X.f , [y.-Y,f 

2 



ln{L) — const — i N 



(A3) 



and the solution effectively involves solving for Xi and 
9j i.e. N + m unknowns, often using a mixture of a direct 
parameter search on ln{L) in conjunction with the following 
differential constraints 



dln{L) 
~dX~ 

dln{L) 
89, 



X,) 



+ 



E 



89, 



dj_ 

8X, 

- f{Xi I 9) 



- f{x. ! 9) 



(A4) 



(A5) 



We illustrate the use of this technique to find the inter- 
cept using the straight line model defined by Yi = aXi + b. 
In this case the differential constraint on Xi (equation IA4|) 
simplifies to 

Xi a{yi — aXi — b) 



Xt 



= 



(A6) 



which implies for a given parameter pair (a, b) that Xi 
is determined by 

ab + XiXj ^^^^ 



X, = 



ayi 



a2 + Xf 



where we have defined A,^ — -P^. 

Given Xi and a trial pair (a, 6), Yi is defined by the 
model and it is then straightforward to compute the log- 
likehood function using a direct search on a suitable grid of 
values (a, b). The result of doing this for the data presented 
in Fig. |X]is shown in the log-likelihood contours in Fig. El 



APPENDIX B: SQL 

An example WFCAM Science Archive SQL query for the 
selection of repeat /f-band observations in the Large Area 
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Survey. Note that the sl.sourcelD < s2.sourceID qual- 
ifier is used to ensure that each match appears only once in 
the output results (and not twice as would happen other- 
wise). 



SELECT sl.kAperMagS as magi , s2 .kAperMag3 as mag2 

FROM lasSource sl.lasSource s2, 

lasMergeLog 11 , lasMergeLog 12, 
lasDetection dl , lasDetection d2, 
Multif rameDetector f 1 ,Multif rameDetector f2, 
CurrentAstrometry cl .Current Astrometry c2 
lasSourceNeighbours x, 

WHERE sl.sourcelD = x.masterObjID AND 
s2.sourceID = x.slaveObjID AND 
si .frameSetlD = 11 . f rameSetlD AND 
s2.fraineSetID = 12.fraineSetID AND 
si .kAperMagSErr < 0.12 AND 
s2 . kAperMagSErr < 0.12 AND 
sl.kclass = -1 AND 
s2.kclass = -1 AND 

fl .multif ramelD = cl .multif ramelD AND 
f2.multiframeID = c2 .multif ramelD AND 
fl.extNum = cl.extNum AND 
f2.extNuin = c2.extNum AND 
dl .multif ramelD = fl .multif ramelD AND 
d2. multif ramelD = f 2 .multif ramelD AND 
dl.extNum = fl.extNum AND 
d2.extNum = f2.extNum AND 
dl.objID = sl.kObjID AND 
d2.objID = s2.k0bjID AND 

11. kmfID = fl. multif ramelD AND 

12. kmfID = f 2. multif ramelD AND 

11. keNum = fl.extNum AND 

12. keNuin = f2.extNuin AND 

si .frameSetID <> s2 . f rameSetID AND 

sl.sourcelD < s2.sourceID AND 

distanceMins < (0.3/60.0) AND 

distanceMins IN ( 

SELECT MIN (distanceMins) 
FROM lasSourceNeighbours 
WHERE masterObjID=x.masterObjID 
) 



where ZP is the zeropoint for the frame (keyword: 
MAGZPT in the FITS header), / is the flux within the 
chosen aperture (e.g. column: APER_FLUX_5), t is the 
exposure time for each combined integration (keyword: 
EXP_TIME), and A is the appropriate aperture correction 
(e.g. keyword: APCOR5). The next term deals with the dis- 
tortion correction caused by the varying pixel scale, where 
(/cor//) is the correction and is derived in Equation [1] The 
final term deals with the extinction correction, where k is 
the extinction coefficient (EXTINCTION) and is equal to 
0.05 magnitudes/airmass in all filters, and x is the airmass 
(keywords: AMSTART, AMEND). 



APPENDIX D: RELEVANT SOURCE AND 
IMAGE PARAMETERS 



APPENDIX C: CONVERSION OF FLUX INTO 
MAGNITUDES 

The processing philosophy is to preserve the image and cat- 
alogue data as counts, and to document all the required 
calibration information in the file headers. Thus recalibra- 
tion of the data requires only changes to the headers, and 
these headers can be reingested into the WSA without the 
need to reingest the full tables. For readers accessing the 
flat files (catalogues and images) rather than the WSA 
database products, we document the methods for converting 
the fluxes into magnitudes and calibrating the photometry. 

m = ZP~ 2.5logio{j) - A - 2.5/ogio(%^) - k{x - 1)(C1) 
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Table Dl. Table of source parameters generated by the WFCAM pipeline and written to the FITS catalogue products, and an accom- 
panying short description for each. 



\ 




T-? 1 1 Tl T» 1 Tl Tn 1 TTl rl/^T" TT^T" 0''ltZO T*rif/^T*OTlr^r^ ITl G"^"T"l/^^" i"\Y'i'\f^y r\T 1 TVl CI (TO ^^"lO'^"0^'■^"li^T^G 
IVLlliillli^ lil_llliU(_-i iUl L.tLOL^ Wi 1 J ill 3LlH_-L Ul t-iL^I (Ji lllidH,"^ (_1.L. LL^l.- L lUilO ■ 


2 


Tc!r»nn<Tl"nl flnv 


^ 1" n VI H 51 TTl rir>fi Tl 1 ("i(~iTi (~it t;iiTYiTTir>ri fliiY wfi1"niTi Ho1"0/^("1(~iti iQ(~iT^lin't"0 tiT^tiT'i" tt"(~itti fio("r>r''t"ir»Ti 

kj L<Jj11L1<X1 LI L1\_.1111IL1*_J11 \J1. oLlllllll\_.Ll 11 LUv WlLlllll *_11_>LL>\_>L1*_J11 loWlJllLJLiL,, C1'IJ(X1 L 1J.W111 Ll\_.L\_.L^LiIl_Jll 

filter is used to cieflue pixel conneetivity cind henee which pixels to include. This helps 
to reduce edge eflects for cill isophoto-lly derived pcircimeters. 


3 


\j\J\JL KA 


TTi+PTicii1"\r— wrpi crn+ori ic;rkT^riri+al r'PTi+T'P— i^t— OTsnT'it'^/' in 

AllliClliSl tj' WClglltCJ. lOL/^ll'J Udl t^CllUl C tjl gl dV lUj' 111 w. 


4 


LLIL L \JL 1X1 -/V 


1— JO 111 111 ell IjC! \JL V/CllliUJHJ Cll IJl • 


5 


I UUOI Ki 


lllttJlloltj' WclglltCU loOpilUtdl t/Clltl fcJ Ol gl dVlLy 111 I . 


6 


Error in Y 


Estimate of centroid error. 


7 


vjxoiiiooiciiii oi^iiicii 


T^pTivprl TrrtTTi thp thrpp intpTTJitv-wpi chtpr) sippOTiH rnornpntsj Tlip pmiivulpTifP to a 

I -/\7-i 1 V v-vl 11 Will Ullv^ Ulll lllLv^llOl v V W ^ICjllU^VJ Ov^^WllLl 1111^111^111* 0» X 11^ Lll VCIilv^lH-'^ XjKJ Ob 

rron ovf^liG^i^n olliT^"^~ir*al l-^^^^^GGl^iTl /HiG'i'T'i tI"i f^n i g 1 1 Go/n r\ ft ot'it^o* 
gclicl dllbtJLl tJiiip lieoii vjdUDDlelll LlloU lULlHUIl lo UocCi tU LltJl 1 Vc. 

GcLUSsio-u sigms, ^ '^b^^^^ 




Ellipticity 


ellipticity ^ 1.0 — (Tdjo'^ 


9 


Position 3,ngl6 


position (iiigle = tingle of ellipse ma.jor 3xis wrt x cixis 


10-16 


ArG&l profile 1 


IVnrnVtPV of ni'ifplci aKovp a ^ptips of thrpshold IpvpI^i vplativp to lopal ciVv TiPvpI*? arp 

L > LllllL^C^l \JL L/lJv^lO <JlVJKJ \ ^ cX OCl 1^1 Ulll v^Ollwl^ 1^ V^IO l^ldljlV^ \jyj IL/^dl iDlv V . 1 V v^lO Ctl C 




XXI ^Ctl k/1 ^111^ ^ 


sipt at T 9T 4T ST 198T whprp T is; thp thrpdholH Thpcjp ran hp thomrht of as; a 
ocij dill J-) ^j-j 'ij-j o J- ... x^o J. w iici c J- lo iiiic iiiii CDinjiu.. J- iicDC i^dii uc (jiii^UKii'-' '-'1 cm d 




Z\ 1 TAyf^Ti 1 o X 
.rVi tJdi JJl UiiifcJ O 


T^r^z-vT* TTl QTi'g T'a/^ial T^Vilfi 1^ ThV^T* /^^hloTl 1 O /^TfOvl Q T^T~\1 Tl CT 1 TTl a fTOG /^Tll"\7 t"ho "fi T'ct aVOH 1 

pUwl llldll o IdiUldl piOlllc. JTOl UcUlCllUtJU, i.e. UVtJl Ictppillg lllldgUo, Ulliy tilt! llloL dltiJdl 




Z\ VO^l 1 T^f rf^Tl 1 O /I 
jT.1 CCtl L/1 LIlllC 


T^T'f^Tl 1 1 G r^rWW TAll t ^ITIM trio T'OGt ^IT'O Got t _ 1 

piUllltJ lo t/WlIip U ttJU dillLl lillfcJ IfcJot dlfcJ oct tU ±. 




A Tpa 1 rirnfi 1 p S 

Jil xUCtX h'-'' ^111^ t-r 






AfPal nTr^Tilp f\ 
JT.1 Cdl Ml MlllC \J 






Arccil profile 1 




17 


Avonl T^T"(~irilr> S 

ii-J. t_ cl 1 LJl WlllS-^ o 


Pov rilpTiHori iTYincro'S tnit; T^piT'nTYir>tr>T' iq hqoH 1"0 flncr trio Qtavt of trir> GPOiiPTir'P of trip 

± Ul LJlt>llV_i.l_^\_l. llll<jjtil_^iz) LrlllO IJ<X1 dllll-^ Lll^l UoL^Li. \j\J lldti Lrlll.^ o Ldl L LJl LrllL^ OCIJ Lldll^C L^l ullC 

deblended components by setting the first in the sequence to 


18 


Pefik height 


In counts relative to local value of sky - also zeroth order aperture flux 


19 


Error in pccik height 




20—45 


A'T^PT'l'nTP flnY 1 

jT-^ClliUlC 11 J- 


T^riPGP ayp a gptipg of HittPi'PTit yaoniG Goft— poctoH aT^pytiifPG opgicttipH to aoponatplv Gam— 

J. llCiSC dl C d OCl ICO \jy UlllCl Cllli 1 dUlUO iSIJlli CUgCU d^JCl li HI Co VJCOlgjllCU \j\J dUCV^UdliClj' iSdlll 




n.TTr»T m flnv 

1 J 1 L \JL 111 11 1 1 A 


T^l o t rio /^n T'"\70— f^T— rrT'f^TTETt" ri r\x t rio tti q if^vit"\7 r\\ i tti a rroG ti m t/~i T^v/^Tfi Ti 'va/H_gi r70/H ^ T^OT't n t'o 
jjlc blltJ t/Ul Vc^Ul gl UW til Ul tlltJ llldiJUllby Ul lllldigco clllU to piUVlUc llXtJU-ol^cU dipcltUlc 




Anpri'iirp fliTy ^ 

L/^l \j \X1. C 11 LIjV ^ 


■niiv-pG for all imao'p^ Thp ^palp ^i^iP foT thp^p ariPTtiiTp^ ciplpptpd hv HpfiniTiP" a ^palp 

11 LIjVV^O 1\^1 dll lllldc^CO. J- lie OVydlC OltiC ll_ll lillCOC dL/Cl U Lll CO lO OCIC*^ UCVl "-'j' ^ClllllllKj d o^^ctic 




1 J 1 1 ^ji 111 11 tiA 


varlint; o^<^P\AAl-rA/r"^ foT *;itp-l-iTi«itTiTmpnt Tn tlip r*a*;p of \AAP(^ A IVT tVii*; " r*OTp" rarliim 

IdWlUO '^^^i VV 111 VI 11^1 OlUCT^lllOUl UllllCllL. Ill UllC ^dOC Ijl VV 1- V^J^IVI LlllO ^wl C Id^lUiO 




Aperture fiux 3 


(rcorc^ has been fixed at 1 . arcsec for convenience in inter- comparison with other 




Error in flux 


datasets. A 1.0 arcsec radius is ecjuivalent to 2.5 pixels for non-interleaved data, 5.0 




Aperture fiux 4 


pixels for 2x2 interleaved data, and 7.5 pixels for 3x3 interleaved data. In ~1 arcsec 




T-l , W/^V* 1 VI "H 1 1 "V 
HjLLKJl 111 ii LIJL 


GO£^l VI Cf QTl y^CrWa VCJ/HlllG Q T^OT*t 111**^ i^rf^TlttllTlG Vrf^l irTlllx/" V/X J^T trlO t /^t Q 1 Tl 1 1 "V C\\ attiilQV 1 TTI — 

occiiig dll 1 c 1 dLii uo duciuuiic ijL'iiijdiiio njui&iiiy tiic uvjudi iilla ui ouciidi iiii 




Aperture fiux 5 


ages. (The rcore parameter is user specifiable and hence is recorded in the output 




Error in flux 


catalogue FITS header.) 




iT-ptJl t HI fcJ 11 UX \j 


The aperture fluxes Eire sky-corrected integrals (summations) with a soft-edge (i.e. 




T^.i'Y'fYt' m flnv 

1 J 1 1 \JL 111 11 l.i-A- 


OTO— Tata flnY onriGion foT oonnoa'nT' oiypIc i r^o^J(7■p\7■p■^ tot o\7PTla'0'r»incr imacrPG trip Aiiypg 

^l-KJ i. OAiOi llli^ UlVlOHJll lljl LHJUllUCLl J ^1 AClO / . ±J.IJWCVClj HJl ^JVCX Id^^Jlllg lllld^CO lillC IIU^CO 




Aperture flux 7 


are derived via simultaneously fitted top-hat functions, to minimise the eflfects of 




Error in flux 


crowding. Images external to the blend are also flagged and not included in the large 




Aperture flux 8 


radius summations. 




H .T'T'J^V 1 "H 1 1 "V 


At^OV'HiIT'O rillV -K 1G 1 V'^r VM Yi TVl OTlriOri IT CITltrlo Tin TYl r^OT* "} C T'O/^n 1 VO/H \' r\ T'OT^T'OGOTlt t ri o fl n "V 

L Lll 11U.J\. O 13 l(_-l„WllllliLliU(_-U Ii tX OlllgiL llLilllUL^l lo ICUUllCU \j\J ICL/lCoCllli tllC llUUl. 




A r^tivi" 1 1 "P^^ tHiiv Q 


for ALL images - this aperture has a radius of rcore. 




w/^f 1 Tl "H 1 1 "V 
Ujl 1 L*! Ill ii LlA. 


^^"QV'^1Tl^T \3^^il"ri T\Q V Q Tn ot */0 trit^ VQ/Hii qvo" fli 1 1 ^ \^ y^crwc^ f*/l 1 /a/VN/t*o/^t*o r-^i vcrw*^ 

kJ Ldl Llllg W 1 Lli JJdl tXlilCbCl ZjVj LllC 1 cLLlll dl C IJ-J \^ ) i- 1 \ ^ i- L-Wi tJ, V '-'/ ^ '^Lll Cj 




Aperture flux 10 


(4) \/2xrcore, (5) 2xrcore, (6) 2\/2xrcore, (7) 4xrcore, (8) 5xrcore, (9) 6xrcore, 




Error in flux 


(10) 7xrcore, (11) Sxrcore, (12) lOxrcore, (13) 12xrcore. 




Aperture flux 11 


Note 4xrcore contains ~99% of PSF flux. 




Error in flux 


The apertures beyond Aperture 7 are for generalised galaxy photometry. 




Aperture flux 12 


Note larger apertures are all corrected for pixels from overlapping neighbouring im- 




Error in flux 


ages. 




Aperture flux 13 


The largest aperture has a radius 12xrcore ie. ~24 arcsec diameter. 




Error in flux 


The aperture fluxes can be combined with later-derived aperture corrections for gen- 
eral purpose photometry and together with parameter 18 (the peak flux) give a simple 
curve-of-growth measurement which forms the basis of the morphological classifica- 
tion scheme. 


46 


Petrosian radius 


Tj, as defined in Yasuda et al. 2001 AJ 112 1104 


47 


Kron radius 


Tk as defined in Bertin and Arnouts 1996 A&A Supp 117 393 


48 


Hall radius 


Th image scale radius eg. Hall & Mackay 1984 MNRAS 210 979 
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49 
50 
51 
52 
53 
54 
55 
56 
57 
58 

59-60 



61 
62 



Petrosian flux 
Error in flux 

Kron flux 
Error in flux 
Hall flux 
Error in flux 
Error bit flag 
Sky level 
Sky rms 
Child/parent 

RA 
DEC 



Classification 
Statistic 



Flux within circular aperture to fe x r^; fe = 2 
Flux within circular aperture to k X r^; k = 2 

Flux within circular aperture to fe X r/d fe = 5; alternative total flux 

Bit pattern listing various processing error flags 
Local interpolated sky level from background tracker 
local estimate of rms in sky level around image 

Flag for parent or part of deblended deconstruct (redundant since only deblended 
images are kept) 

RA and Dec explicitly put in columns for overlay programs that cannot, in general, 
understand astrometric solution coefficients - note r*4 storage precision accurate only 
to 50iuas. Astrometry can be derived more precisely from WCS in header and XY in 

parameters 5 & 6 

Flag indicating most probable morphological classification: eg. -1 stellaj:, +1 non- 
stellar, noise, -2 borderline stellar, -9 saturated 

An equivalent N(0,1) measure of how stellar-like an image is, used in deriving pa- 
rameter 61 in a "necessary but not sufficient" sense. Derived mainly from the curve- 
of-growth of flux using the well-defined stellar locus as a function of magnitude as a 
benchmark (see Irwin et al. 1994 SPIE 5493 411 for more details). 



Table D2. Relevant photometric parameters measured by the pipeline and written to the FITS headers. These values are computed 
per-detector and stored in the headers for each image and catalogue extension. The names by which the parameters are stored in the 
WFCAM Scionre Arrhi\'e tables are also gi\'cii. 



FITS keyword 


WSA parameter 


Description 


AMSTART 


amStart 


Airmass at start of observation 


AMEND 


amEnd 


Airmass at end of observation 


PIXLSIZE 


pixelScale 


[arcsec] Pixel size 


SKYLEVEL 


sky Level 


[counts/pixel] Median sky brightness 


SKYNOISE 


skyNoise 


[counts] Pixel uoisc at sky level 


THRESHOL 


thresholdlsoph 


[counts] Isophotal analysis threshold 


RCORE 


coreRadius 


[pixels] Core radius for default profile fit 


SEEING 


seeing 


[pixels] Average stellar source FWHM 


ELLIPTIC 


avStellarEll 


Average stellar ellipticity (1-b/a) 


APCORPK 


aperCorPcak 


[magnitudes] Stellar aperture correction - peak height 


APCORl 


aperCorl 


[magnitudes] Stellar aperture correction - core/2 flux 


APCOR2 


aperCor2 


[magnitudes] Stellar aperture correction - core/\/2 flux 


APCOR3 


aperCor3 


[magnitudes] Stellar aperture correction - core flux 


APCOR4 


aperCor4 


[magnitudes] Stellar aperture correction - \/2x core flux 


APCOR5 


aperCorS 


[magnitudes] Stellar aperture correction - 2x core flux 


APCOR6 


aperCor6 


[magnitudes] Stellar aperture correction - 2^/2x core flux 


APCOR7 


aperCor7 


[magnitudes] Stellar aperture correction - 4x core flux 


MAGZPT 


photZPExt 


[magnitudes] Photometric ZP for default extinction 


MAGZRR 


photZPErrExt 


[magnitudes] Photometric ZP error 


EXTINCT 


extinctionExt 


[magnitudes] Extinction coefliicient 


NUMZPT 


numZPCat 


Number of 2MASS standards used 


NIGHTZPT 


nightZPCat 


[magnitudes] Average photometric ZP for the fllter for the night 


NIGHTZRR 


nightZPErrCat 


[magnitudes] Photometric ZP a for the fllter for the night 


NIGHTNUM 


nightZPNum 


Number of ZPs measured for the fllter for the night 



